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^Bounding  theorems  for  the  changes  in  elastic  energy  due  to  a  crack  in  a 
finite  body  and  crack  tip  Stress  Intensity  Factors  are  established  by  use  of  the 
variational  principles  of  the  theory  of  elasticity.  It  is  shown  that  the  change 
in  elastic  energy  due  to  a  crack  in  a  finite  body  is  larger  than  the  change  in 
energy  due  to  the  same  crack  in  an  infinite  body  with  the  same  boundary  condit¬ 
ions  when  stresses  are  prescribed  on  the  external  boundary,  and  smaller  when 
displacements  are  prescribed.  The  energy  changes  can  be  expressed  as  functions 
of  the  crack  tip  Stress  Intensity  Factors  and  for  special  loadings.  (over) 


DD  t  1473  EDITION  OF  I  NOV  SS  IS  OBSOLETE 


SECURITY  CL  ASSIFICA'ION  OF  This  PACE  fl*»i»/l  Dai«  Ento  rod) 


Abstract,  continued 


“bounds  for  single  mode  Stress  Intensity  Factors  are  obtained.  The  obtained 
inequalities  are  in  agreement  with  known  numerical  solutions  of  finite  cracked 
bodies. 
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1.  Introduction 

The  present  work  is  concerned  with  new  bounding  theorems  in  fracture 
mechanics.  The  theorems  are  established  by  the  use  of  the  variational 
Principles  of  the  Theory  of  Elasticity.  The  theorems  show  the  influence  of 
the  finite  size  of  the  cracked  body  on  the  change  in  elastic  energy  due  to 
the  crack  as  a  function  of  the  prescribed  boundary  conditions.  For  special 
loadings,  when  single  inodes  occur,  it  is  shown  that  the  stress  intensity 
factors  (abbreviated  SIF)  of  a  finite  cracked  body  with  prescribed  stresses 
on  the  external  boundary  will  always  be  larger  than  the  crack  tip  SIF  of  am 
infinite  cracked  body  with  the  3ame  boundary  conditions  and  that  the  SIF  of 
a  finite  cracked  body  with  prescribed  displacements  on  the  external  boundary 
will  be  always  smaller  than  the  SIF  of  am  infinite  cracked  body  with  the 
same  boundary  conditions. 

2.  The  Bounding  Procedure 


The  problem  of  bounding  the  change  in  elastic  energy  due  to  a  crack  in 
a  finite  body  will  here  be  considered  in  terms  of  the  classical  extremum 
principles  of  minimum  potential  and  complementary  energies.  These  principles 


have  been  previously  used  for  bounding  the  effective  elastic  moduli  of 
heterogeneous  materials  [l]  and  recently  for  bounding  effective  moduli  of 
cracked  material  [2]  . 


t'  ,  ■'  •;  ,i 


ft 


OO  \ 
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Finite  BcflY  with  prescribed  Stresses. 


The  change  in  energy  due  to  a  crack  in  a  finite  body  with  prescribed 
stresses  on  the  external  boundary  will  be  bounded  by  use  of  the  principle 
of  minimum  potential  energy.  Consider  the  finite  body  of  volume  V,  con¬ 
taining  a  stress  free  crack  (fig.  la).  Traction  boundary  conditions  are 
prescribed  on  the  external  boundary  . 

Ti(S1)  =  T*  (1) 


The  potential  energy  of  the  cracked  body  is 


UP1  "  2 


ijl  tijl 


dV  - 


li  il 


dS 


(2) 


where  and  u ^  are  the  unknown  stresses,  strains  and  displacements 

in  the  body.  The  potential  energy  stored  in  the  cracked  body  is  given  by, 


PI 


“p,  -  ®, 


(3) 


where  is  the  potential  energy  of  the  body  of  volune  V1  subject  to  (1) 
in  the  absence  of  the  crack,  and  is  a  function  of  the  stresses 
and  strains  in  the  uncraoked  body, 


dV 


(4) 
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and  SU  is  the  change  in  potential  energy  due  to  the  crack  and  according  to 
a  well  known  result,  Eshelby  C 3 ] ,  can  be  expressed  as 

«u  .  i  aij  [  <5) 

on  the 
crack 

where  [u^]  is  the  displacement  jump  accross  the  crack  faces.  From  (5)  it 
is  seen  that  only  the  displacement  in  the  immediate  vicinity  of  the  crack 
has  to  be  known  to  compute  <SU. 


The  potential  energy  functional  is  defined  by: 


°P1  *  2 


aij  £ij  dV  ‘ 

V1 


T l  ut  dS 


(6) 


where  u^  is  an  admissible  displacement  field  which  must  merely  be  continu¬ 
ous  in  the  region  excluding  the  crack.  Since  in  the  present  problem  boun¬ 
dary  tractions  are  prescribed  u^  is  not  restricted  by  boundary  conditions. 
The  "stresses"  and  "strains",  cr^  and  are  defined  by 


*it  !i  (ui.j  *uj.u  <a) 


(7) 


aij  *Cijkl  ®ia 


(b) 


The  extremum  principle  states  that 
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0P1  *  UP1  (8) 

Thus  the  problem  of  bounding  reduces  to  that  of  constructing  an  appro¬ 
priate  admissible  displacement  field.  To  do  this  consider  another  larger 
finite  body  of  volume  V  including  and  containing  the  same  crack,  where 


V  s  V1  +  V2  >  V,  (9) 

a 

and  apply  the  tractions  to  the  external  boundary  S2 

W  s  *i  (10) 

The  displacement  solution  to  this  problem  is  The  part  of  this  dis¬ 

placement  field  contained  in  the  region  of  the  larger  body  is  an  admis¬ 
sible  displacement  field  for  the  body  of  volume  subjected  to  (1),  since 
it  is  continuous  within  the  volume  V^,  excluding  the  crack. 

=  ui2(in  ?.,)  (11) 


The  displacement  field  u^2  can  be  expressed  as  the  sum  of  the  displacement 
u°2  in  the  body  V  without  the  crack  and  the  perturbation  u|2  due  to  the 
presence  of  the  crack. 


Ui2  * 


12 


i2 


(12) 


The  strains  and  stresses  associated  with  ui2  are 


eij2  3  eij2  +  £ij2 


(a) 


aiJ2  S  0ij2  +  aij2 


(b) 


Thus 


eij  *  eij2(in  V 


(a) 


aiJ  =  aij2(in  V1} 


(b) 


The  potential  energy  of  the  larger  body  is  defined  as 


JP2  2  J  aij2  eij2  dV  "  J  ^i  “i2 

12  2 


T?  dS 


and  can  be  expressed  as 


°P2--2 


Jv+y  ^2^2d7-5D2 

V1  v2 


(13) 


(14) 


(15) 


(16) 


In  order  to  bound  the  real  potential  energy  of  the  finite  body  of  volume  V. 


the  expression  0^  is  defined  by 


®p--H 


e°  .  dV  -  <5u,. 


i12  ij2 


(17) 
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Comparing  (16)  and  (17)  it  is  easy  to  see  that  Op  takes  the  form 


UP  *  UP2  +  2 


°ij2  £ij2  dV 


(18) 


By  substituting  (11)  and  (14)  into  (6)  the  potential  energy  functional  for 
takes  the  form: 


"  2  L  aij2  eij2 
V1 


dV  - 


Ti  ui2  dS 


(19) 


It  will  now  be  shown  that 


°p  sp 


(20) 


from  which  it  then  follows  by  (8)  that 


0  2.  U  , 

P  Pi 


(21) 


Assuming  (20)  to  be  correct  it  can  be  written  in  the  form 


H 


7Jv+v°1J2  '«2 

V1  2 


117  ’  |s  Ti  “12  ^  +  2  j„  °1J2  s«2d,iff. 


eiJ2  dVi7jv,JiJ2elj2dV- 
2  1 

(22) 


Ti  u12  43 


or  after  rearranging 


m|h* 
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(aij2£ij2'H7 


o  o 
ij2£ij2 


)  dV- 


12 


dS+ 


u 


12 


dS  >  0 


(23) 


It  will  now  be  assumed  that  has  the  special  form 


Ta  o 

1  *  ij  “j 


(24) 


both  on  S1  and  on  3^.  This  is  known  as  homogeneous  traction  boundary  con¬ 
ditions  and  the  elasticity  stress  solutions  in  this  case  are 


o  o  o 

aij1  s  aiJ2  s  aij 


(25) 


for  homogeneous  elastic  bodies  of  arbitrary  shape.  The  associate  strain 
and  displacements  are 


o  o  o 

eij1  s  eij2  s  eiJ 


u11  (x)  a  ui2  (x) 


a  S 


ijkl  aij 


(26) 


uiU)  ■  elj  xi 


It  now  follows  from  virtual  work  that 


(  a  ( 

f 

*1  *  1: 

s  TI  Ui2dS  s  -  j 

dV 


(27) 


Introducing  (27)  together  with  03)  and  (25)  into  (22)  we  obtain 
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<  -  T  °li  s(j2  +  i  °'l}2  +  2  ah 2  e;J2>  dV  i 0 


(28) 


By  the  symmetry  of  the  with  respect  to  ij,kl  interchange 


cr3  g»  -  j«  c° 

i4  iJ2  iJ2  ij 


(29) 


and  (28)  reduces  to 


2  aij2  eij2  dV  ^  0 


(30) 


Evidently  (30)  is  correct  as  the  integrand  is  positive  definite  since 
it  is  an  elastic  strain  energy  density.  This  proves  the  inequality  (20), 
therefore  also  the  inequality  (21).  Expressing  0p  by  (17)  and  Upt  by  (3) 
and  (4)  and  using  (25),  (21)  becomes: 


dV 


Su2  z- 


if  a0 


and  therefore 


(3D 


6U1  i  (32) 

in  words:  gdaa  Ida  sans,  daaggaaaaua  .sfcraaa  Eaundacy  conditions  ac a 
ar.gacridgd  an  jma  asternal  teuadary  a£  iaa  finite  Sadias.  Ida  aaam  change 
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du&  la  Ha.  flcaacaas.  zL  a  crack  will  i&  larger  la  Ha  saallac  H&u 

A  limiting  case  of  the  problem  considered  is  when  the  larger  body  of 
volume  V  becomes  infinite.  It  follows  that  the  change  in  energy  due  to  a 
crack  in  a  finite  body  subjected  to  stress  boundary  conditions  is  always 
larger  than  the  energy  change  is  an  infinite  body  subjected  to  the  same 
boundary  conditions.  Thus  Eq.  (32)  can  be  extended  to  read 

«01  1  <SU2  1  <5U°°  (33) 

Finite  body,  with  prescribed  displacements. 

The  change  in  elastic  energy  due  to  a  crack  in  a  finite  body  with 
prescribed  displacements  on  its  external  boundary  will  be  bounded  by  use  of 
the  principle  of  minimum  complementary  energy.  The  procedure  is  very  simi¬ 
lar  to  that  followed  above,  but  instead  of  am  admissible  displacement  field 
an  admissible  stress  field  has  to  be  constructed. 


Let  us  prescribe  homogeneous  displacements  on  the  external  boundaries 
S1  and  Sg,  of  the  two  finite  cracked  bodies  of  volumes  V1  and  V=V1+V2  il¬ 
lustrated  in  fig.1. 


ut<V  5  eIj  XJ 


(a) 


(34) 


u?(S2) 


'ij  J 


Cb) 


•  v 
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The  stress,  strain  and  displacement  fields  in  the  two  bodies  are 

aijr  un  and  aij2  £ij2’  ui2'respeotively  the  Presence  of  the 

crack  and  a°j1«  eiji»  and  aij2*  eij2’  ui2  in  the  absenee  of  the 
crack.  The  complementary  energies  of  the  two  bodies  are 

UC1  s  2  l  aijl£ijl  ^  '  [  Til  ui  dS  "  “  2 
V1  s 

UC2  *  2  L  _  aij2£ij2  ^  “  L  Ti2  Ui  ^  "  ”  I  L  aij2eij2  dV  +  6U2  (36) 

V1  v2  a2  vi 

where 


ai11£ijl  dV  + 


«U. 


(35) 


Tis 


aij  "j 


(37) 


As  in  the  previous  case  the  energy  changes  6U  are  given  by  (5). 


The  stress  solution  to  the  larger  body,  of  volime  V,  over  the  part  V1 
of  this  volume  is  an  admissible  stress  field  for  the  problem  of  the  body 
V^,  since  it  satisfies  equilibrium  everywhere  in  and  has  no  displacement 
boundary  conditions  to  satisfy,  thus 


°i32  (1“  V 


(38) 
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The  potential  energy  functional  for  the  body  V1  is  defined  as 


°C1  "  2  „  aij  gij  dV 


Ti  ui  48 


(39) 


where 


T.  =  a..  n. 

1  ij  J 


(40) 


and  is  related  to  a  by  (7b) . 


The  principle  of  minimum  complementary  energy  states  that 


UC1  *  UC1 


(41) 


Define  0C  as 


■  -  2  L 


°«2  eIj2  dT  +  SD2 


(42) 


From  (36)  it  follows  that  D.can  be  expressed  as 


Sc",C2+l{. 


aij2eij2  dV 

v2 


H 


Vv2°1J2SlJ2  dV  'S2 


-  L>  < 


dS  + 


(43) 


.0  -O 


+ 1  (v  ^:3i « 

v2 


Page  14 


For  boundary  conditions  (34),  (25-26)  are  true  for  elastic  homogeneous 
bodies  of  arbitrary  shape.  Expressing  the  admissible  stresses  and  the  re¬ 
lated  strains  by  (13)  and  using  the  principle  of  virtual  work  for  the  vo¬ 
lume  it  can  be  shown,  as  in  the  previous  case  that 

°c  i  “c  <«> 

It  then  follows  that 

5Ci0C1  (45) 

From  (35)  and  (43)  substituted  into  (45)  we  obtain 

6U2  1  60,  (46) 

In  words:  jflaga  lUfi.  same  fagflOggagPUg  displacement  boundary  conditions  gee. 
araagrlbed  aa  iaa  asternal  bauadarlea  of  two  finite  bodies,  the  energy 

change  due  to  the  presence  of  g  crack  will  be  smaller  in  the  small  an  body. 

Again,  a  limiting  case  is  obtained  when  V  becomes  infinite:  the 
change  in  elastic  energy  due  to  a  crack  contained  in  a  finite  elastic  body 
with  prescribed  displacements  on  its  external  boundary  is  always  smaller 
than  the  change  in  elastic  energy  due  to  the  same  crack  in  an  infinite  body 
with  the  same  prescribed  displacements  at  infinity.  Eq.  (46)  extends  to 


read: 
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5U°°  1  6U2  1  60 

3.  Bounding  zL  SUL 


(47) 


In  the  preceding  general  bounds  on  the  changes  in  elastic  en¬ 
ergies  due  to  the  presence  of  a  crack  in  a  finite  body  were  obtained. 
These  bounds  can  be  transformed  into  bounds  on  SIF  for  special  loadings. 


The  change  in  energy  due  to  a  Griffith  crack  imbedded  in  an  isotropic 
body  as  given  by  (5)  can  be  expressed  as  a  function  of  the  crack  tip  SIF. 
For  a  3-dimensional  body 


a 

[Kj(a)  ♦  KjjU)  (1-v)Kjn(a)]  da 
o 


(48) 


When  the  cracked  body  is  2-  dimensional ,  only  Modes  I  and  II  exist  and 
a 

6U  a  jjr-  j  [Kj(a)  ♦  K^(a)]  da  (49) 

o 

where  a  is  half  the  crack  length. 

When  the  cracked  body  becomes  infinite  the  stress  intensity  factors  are 


_  1/2 
KI  a  a22  a 


(a) 


K 


II 


a  a 


12 


a 


1/2 


(b) 


(50) 


*111 


_o  1/2 
ff23  a 


(c) 
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where  cr®^  are  the  tensile  and  shear  stresses  applied  at  infin 

ity,  when  stress  boundary  conditions  are  prescribed,  or 


(a) 

(b)  (51) 

(c) 


when  displacements  of  the  kind  u^0-  £°jXj  are  applied  on  the  external  boun¬ 


dary. 


For  the  general  case  of  mixed  mode  loading,  only  the  combination 
2  2  2 

[Kj  ♦  Kjj  +  (1—  v)KIII]  of  the  stress  intensity  factors  can  be 
bounded.  For  special  single  mode  loadings  bounds  on  single  S1F  can  be  ob¬ 
tained.  A  few  examples  will  be  considered  in  the  following. 


Example  1 : 


a)  Stress  boundary  conditions. 


Consider  a  finite  plane  rectangle  of  dimentions  2b  and  2c  containing  a 
crack  of  length  2a.,  (fig. 2).  with  prescribed  boundary  conditions  on  the 
external  edges  of  the  form: 
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a22(x2  s  *■  c)  =  a° 

an  (x1  s  ±  b)  s  0  (52) 

ai2  (x1  s  *  b)  *  ai2  (x2  5  *  c)  8  0 


Since  the  loading  is  symmetric  this  is  a  mode  I  loading,  and  the  only  SIF 
existent  is  Kj.  The  change  in  elastic  energy  due  to  the  crack  is 

a 

<5U  -  —  |  K^(a)  da  (53) 

o 


The  change  in  elastic  energy  due  to  a  crack  in  an  infinite 
the  same  boundary  conditions  is 

/**?«*• 

o  o 


body  with 


(54) 


Substituting  (53)  and  (54)  into  the  inequality  (33)  we  obtain 

a  a 

(  K*  (a)  da  >  |  K^2  (a)  da  (55) 

o  o 

This  inequality  will  hold  for  any  crack  length  only  when  the  same  inequali 
ty  holds  for  the  integrands,  therefore 


Ki  *  lT 


(56) 
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or 


Kil 


o  1/2 
a  a 


(57) 


b.  Blanlafigaant.  boundary  conditions 


Consider  again  the  rectangle  of  Fig  2.  with  prescribed  displacement 
boundary  conditions  on  the  external  edges  of  the  rectangle  of  the  form: 

u2  (x2  8  *  c)  s  *■  vo 

u-l  (x2  S  ±  d)  =0  (58) 

«2  (*1  s  4.  b)  S  u1  (xt  S  b)  S  0 


This  deformation  is  symmetric  so  only  the  mode  I  SIF  will  enter  the 

energy  change  expression.  Me  can  express  the  prescribed  displacement  v  as 

o 


o 

vo  se  c 


The  stress  intensity  factor  for  an  infinite  body  with  prescribed  displace¬ 
ments  at  infinity  according  to  (51a),  is 


s"  .  Ee°  .1/2 


(59) 
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When  the  changes  in  elastic  energy  for  the  finite  and  infinite  rectangle 
are  expressed  as  a  function  of  the  mode  I  SIF  and  inequality  (47)  is  used 
the  following  bound  on  K^(a)  is  obtained: 


KjU)  i  K*  (a) 


(60) 


or 


Kj(a)  i. 


Ev 


1/2 


(61) 


Isida  [7]  calculated,  numerically,  the  mode  I  SIF  for  a  finite  rectan¬ 
gle  with  various  boundary  conditions.  For  boundary  conditions  (51)  his  re¬ 
sults  satisfy  (57),  i.e.  the  SIF  are  always  larger  than  the  one  for  an  in¬ 
finite  cracked  body. 


A  second  problem  considered  by  Isida  was  that  of  the  finite  rectangle 
with  the  following  boundary  conditions 


u2  <*2  s  *  c)  s  ±  vo  (a) 

Uj  (*2  *  ±  c)  *  0  (b) 

a11  (*•,  *  ±  b)  *  a12  (x1  =  ±.  b)  a  0  (0) 


(62) 


+ 

te. 
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It  should  be  pointed  out  that  the  results  for  this  case  cannot  be 

bounded  by  the  theorems  developed  in  the  present  work.  Since  the  boundary 

conditions  are  mixed  [  (62c)  are  traction  boundary  conditions].  However 

Evo  1/2 

these  results  also  satisfy  (61),  i.e.  is  always  smaller  than  — —  a.  . 

It  should  be  emphasized  that  there  are  no  similar  inequalities  for  the 
SIF  of  a  finite  body  with  mixed  boundary  conditions.  It  may  be  either 
larger  or  smaller  than  the  SIF  for  the  corresponding  infinite  body.  Isida 
[4]  also  considered  the  case  of  the  rectangle  with  the  following  loading 

u2  (x2  a  ±  c)  *  ±  vQ 

C12  (x1  =  b)  3  ai2(x2  =  o)  s  0  (63) 

°ti  (*, » *  0 


For  this  case  the  results  obtained  where  sometimes  larger  than  and 
sometimes  smaller,  depending  on  the  ratios  a/b  and  c/b. 

E-ramnl  m  2 ! 

A  second  example  considered  will  be  that  of  a  finite  rectangle  of  di- 
mentions  b  and  2h  containing  a  cracx  of  length  a,  (fig. 3). 

The  corresponding  infinite  body  for  this  problem  will  be  that  of  a 


i 
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j 

half  plane  containing  a  crack.  The  mode  I  SIF  for  the  infinite  body  is 

K~  =  1.72  j°  a1/2  (64) 

o 

when  a  stress  a  perpendicular  to  the  crack  is  prescribed  at  infinity  and 
ICj  s  1.12  E  £°  a1/2  (65) 

when  a  displacement  field 

u2°  =  e°  x2  (66) 

is  prescribed  at  infinity. 

Various  solutions  to  the  problem  of  the  finite  rectangle  with  boundary 
conditions 

°22  ^X2  *  *■  s  a° 

a12  (*1  S  0,b)  S  a12  (x2  s  t  h)  *  o  (67) 

"it  <*i  *  o-w  1 0 

have  been  summarized  in  [8]  and  the  K^(a)  obtained  are  always  larger  than 
1.12  a°  a1/2  as  predicted  by  eq.  (56), 
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Also  the  mode  III  loading  with  prescribed  stresses  for  a  finite  bar 
with  an  edge  crack  was  solved  numerically,  [5],  and  the  results  are  also  as 

OO 

expected:  the  obtained  (a)  are  larger  than  . 

4.  fianeiuaiaaa 

By  using  the  variational  principles  of  the  theory  of  elasticity  it  has 
been  shown  that  the  change  in  energy  due  to  a  crack  in  a  finite  body  incre¬ 
ases  as  the  dimensions  of  the  body  decrease  and  is  always  larger  than  the 
change  in  energy  due  to  a  crack  in  an  infinite  body,  when  stress  boundary 
conditions  are  prescribed.  When  displacement  boundary  conditions  are  pres¬ 
cribed  the  inverse  is  true.  As  the  dimensions  of  the  body  decrease,  the 
change  in  elastic  energy  also  decreases  and  the  change  in  energy  due  to  a 
crack  in  an  infinite  body  is  always  larger  than  in  a  finite  one. 

The  changes  in  elastic  energy  can  be  expressed  as  a  function  of  the 
crack  tip  SIF  and  for  special  loadings  the  SIF  can  be  shown  to  be  always 
larger  (for  prescribed  stresses)  or  Smaller  (for  prescribed  displacements) 
than  the  SIF  of  a  crack  in  an  infinite  body. 

The  changes  in  elastic  energy  can  be  related  to  the  crack  criticality 
criterion.  The  decrease  in  the  size  of  a  finite  body  will  precipitate  fai¬ 
lure  when  stresses  are  applied  and  postpone  it  when  displacements  are  ap¬ 
plied.  When  alternating  stresses  are  applied  to  a  finite  fatigue  specimen 
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and  a  crack  is  growing  during  the  fatigue  process,  the  crack  tip  SIF  incre¬ 
ases  with  the  dimensions  of  the  crack  and  so  does  the  crack  criticality. 

Existing  solutions  to  crack  problems  in  finite  bodies  have  been  shown 
to  be  in  agreement  with  the  theorems  developed  in  the  present  work. 
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flours  fiacliana 

Fig.  1:  Finite  cracked  bodies  of  arbitrary  size 

Fig.  2:  Finite  rectangle  with  central  crack 

Fig.  3:  Finite  rectangle  with  edge  crack 


